Lattice effects and entropy release at the low-temperature phase transition in the 
spin-liquid candidate k-(BEDT-TTF) 2 Cu 2 (CN) 3 
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The spin-liquid candidate k-(BEDT-TTF)2Cu2(CN)3 has been studied by measuring the uniax- 
ial expansion coefficients on, the specific heat, and magnetic susceptibility. Special emphasis was 
placed on the mysterious anomaly around 6 K - a potential spin-liquid instability. Distinct and 
strongly anisotropic lattice effects have been observed at 6K, clearly identifying this feature as a 
second-order phase transition. Owing to the large anomalies in on, the application of Griineisen 
scaling has enabled us to determine the corresponding specific heat contribution and the entropy 
release. Comparison of the latter with available spin models suggests that spin degrees of freedom 
alone cannot account for the phase transition. Scenarios, involving charge degrees of freedom, are 
discussed. 

PACS numbers: 71.30.+h, 75. 50. Mm, 71.27.+a, 74.70.Kn 
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Organic charge-transfer salts show a rich variety of 
quantum phases resulting from the interplay of elec- 
tronic correlations, low dimensionality and frustrated 
magnetic interactions (see refs. [l[ and Q for recent re- 
views). Recently, much excitement has been generated 
by the proposal of a quantum spin-liquid (QSL) state 
in the k-(BEDT-TTF) 2 Cu 2 (CN) 3 salt [3]. This sys- 
tem is a half-filled Mott insulator, where (BEDT-TTF)J 
dimers form a quasi-two dimensional (quasi-2D) trian- 
gular lattice characterized by mier-dimer hopping am- 
plitudes t' ~ t. The material lacks long-range magnetic 
order down to 32 mK ^3j , a small fraction of the estimated 
nearest-neighbor Heisenberg exchange coupling J/ks — 
250 K [3], [J This observation has stimulated numerous 
works aiming at an identification of the QSL state and its 
potential instabilities. A number of important questions 
has been raised such as the nature of the low-energy ex- 
citations |i[ , as well as the actual degree of frustration. 
According to refs. and H, the system is less frustrated 
(t'/t ~ 0.8) than was previously thought [9(, raising the 
question of how, under these conditions, a QSL state can 
be stabilized over magnetic order 0, llbHl5j . 

A central issue, which may hold the key for under- 
standing the puzzling low-T properties, is posed by the 
mysterious anomaly around 6 K, about which very little is 
known. It appears as a hump in the specific heat the 
NMR relaxation rate [3], and the thermal conductiv- 
ity [||. Various scenarios have been suggested including 
a crossover from a thermally to a quantum disordered 
state [j| or an instability of the QSL. For the latter, a 
spin-chirality ordering 17j. a & vortex formation [18| , 
a pairing of spinons [19 1, a n exciton condensate [20j, or 
a nematic paired state 2l( have been discussed. Unfor- 
tunately, neutron scattering - the technique of choice for 
identifying a potential order parameter - is very difficult 
to apply due to the organic nature of the material. How- 
ever, according to refs. [l9| and [2(j, the lattice strain 



could constitute such a sensitive probe as the QSL insta- 
bilities discussed there break the lattice symmetry. 

In this Letter, we report on a study of lattice effects 
by means of thermal expansion measurements. Our data 
reveal a distinct and strongly anisotropic anomaly at 6 K 
- clear evidence for a second-order phase transition. The 
application of a thermodynamic analysis enabled us to 
determine the entropy release at the transition. This 
quantity provides a crucial test for any model attempting 
to describe the low-T properties of this material. 

The thermal expansion measurements were carried out 
by using an ultrahigh-resolution capacitive dilatometer 
(built after [22J) enabling the detection of length changes 
Al > 1CP 2 A. The calorimetric and magnetic measure- 
ments were performed employing a micro-calorimeter [23[ 
and a SQUID magnetometer (Quantum Design MPMS), 
respectively. The single crystals used were prepared by 
following the standard procedure [lij]. The crystals, 
taken from two independently prepared batches, have the 
shape of flat distorted hexagons with dimensions of about 
1 x 0.8 x 0.1 mm 3 . The uniaxial pressure exerted on the 
crystal by the dilatometer ranges from 1 to 6 bar. 

In Fig.[TJ we show the results of the uniaxial expan- 
sion coefficients a^T) = l^dh/dT measured along the 
in-plane i — b and c axes and the cross-plane a axis be- 
low 200 K. Besides the distinct anomalies at low temper- 
atures, which will be discussed below, the data reveal 
strongly anisotropic expansivities. Surprisingly, a partic- 
ularly strongly pronounced anisotropy is found for the 
in-plane data and a Cl both of which deviate markedly 
from an ordinary lattice expansion. Upon cooling, a c be- 
comes reduced almost linearly down to 70 K, changes sign 
at 50 K and passes through a shallow minimum around 
30 K. This contrasts with the pronounced maximum in 
at at T max ~ 70 K. The temperature range around T max 
is also distinct in the out-of-plane data a a , where a shoul- 
der appears. The anomalous T-dcpcndencics in at sug- 
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FIG. 1: (Color online) Uniaxial expansivities an of k-(BEDT- 
TTF)2Cu2(CN)3 along the in-plane i = b and c axes (left 
scale) and along the out-of-plane i — a axis (right scale). 
Inset shows the volume expansion coefficient ft. 



gest that, besides phonons, other excitations contribute 
significantly to the thermal expansion. For example, the 
sign change in a c for T below about 50 K followed by 
a minimum indicates a substantial negative contribution 
which is maximum around 30 K. Such a broad negative 
anomaly, lacking any signature in the magnetic prop- 
erties [3j, may have different origins such as geometri- 
cal frustration and/or quenched disorder (see, e.g. [H]). 
Likewise, a rounded feature, similar to the one observed 
in the magnetic susceptibility x{T) at T x w 85 K Q, 
is expected at T a ~ T x as a result of short-range anti- 
ferromagnetic correlations [26|. Irrespective of the na- 
ture of the various anomalies in a,, the distinct ab vs. a c 
anisotropy implies pronounced T-dependent in-plane lat- 
tice distortions. The effect is particularly distinct for 
T < 50 K, where upon cooling the &-axis lattice parame- 
ter strongly contracts (large positive ab) while the c-axis 
lattice constant expands (a c < 0). Since the hopping am- 
plitudes t' and t depend sensitively on the lattice param- 
eters, we expect that cooling in this temperature range 
is accompanied by an increase of t'/t (cf. inset of Fig. [5]). 

Turning to the anomaly at 6K, shown in Fig. [2] on ex- 
panded scales, the uniaxial expansivities reveal a distinct 
peak of positive (a c ) and negative (ab and a a , cf. Fig.QJ 
sign, which is most strongly pronounced in a c . The shape 
of this feature and its sharpness are clear indications of a 
phase transition, albeit of distinctly non-mean-field type 
[27j . signalizin g th e presence of strong critical fluctua- 
tions (see, e.g. 28J, |29|). Measurements taken upon cool- 
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FIG. 2: (Color online) In-plane expansivities of k-(BEDT- 
TTF)2Cu2(CN)3 on expanded scales with a c taken in B = 
(red circles) and 8T (green circles). Inset: 2D triangular- 
lattice dimer model with hopping amplitudes t' and t. 



ing and warming at a very slow rate of ± 1.5 K/h failed 
to detect any hysteresis, consistent with a second-order 
transition. Around 3K the a c data reveal indications for 
yet another anomaly of much smaller size, reproducible in 
detail in consecutive runs. We stress that the features at 
6 K and 3 K in a c remain unaffected by a magnetic field 
of 8T applied along the c-axis (cf. Fig. [2]). Remarkably 
enough, the volume expansion coefficient /3 = a a +at,+a c 
shows a much less peculiar behavior (cf . inset of Fig.QJ . 
In particular, /3 varies smoothly around 50 K and lacks 
any anomaly at 6K, i.e., A/3 \qk~ 0. According to the 
Ehrcnfest relation, this implies that the 6K transition is 
practically insensitive to hydrostatic pressure. 

In Fig. [3] we show results of the specific heat on a 
small single crystal of mass < 100 fig for temperatures 
2K < T < 10 K. The data reveal a smooth increase with 
T and a hump-like feature around 6 K, consistent with 
literature results Evidently, the quantity of interest 
- the contribution associated with the phase transition 
SC trans - is difficult to separate from the background C bg . 
Although C hg is likely to be dominated by the lattice spe- 
cific heat C ph (unfortunately unknown) it may also con- 
tain other contributions. Attempts to separate SC trans 
by subtracting from the measured specific heat C ph of 
a related salt with a different anion [f|, involve consid- 
erable uncertainties. In order to overcome this problem, 
we use an Ansatz, which has proved particularly valu- 
able for analyzing phase transitions in organic materials 
[30j ] . The approach is based on the assumption of a pro- 
portionality between corresponding contributions to the 
expansivity and specific heat Sai oc SC. This so-called 
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FIG. 3: (Color online) Specific heat of k-(BEDT- 
TTF) 2 Cu 2 (CN) 3 between 2 and 10 K. Solid line represents 
the result of a least-squares fit of Cfu(T) (see text), yielding 
a smooth background C 9 (broken line). Inset: entropy S arl 
obtained by integrating 5C an /T from 2K to 10 K. 



Griineisen-scaling applies provided that there is only one 
relevant energy scale in the T range of interest (see j2(| 
and references cited therein). In what follows, we focus 
on the range 2K < T < 10K, where this assumption 
should be valid. First, the data in Fig. [2] were used to 
extract the anomalous contribution Sa^ n associated with 
the peak anomaly around 6K. To this end, smooth back- 
ground curves were subtracted from a c and ctb such that 
the remaining contributions show the scaling behavior 
5a^ n oc Sa^ n required by thermodynamics (cf. Fig.[4j. 
Note that due to the large size of the anomaly in a c , 
uncertainties in the background have only little effect on 
the resulting 8a® n . With the so-derived 5a°: n a least- 
squares fit to the specific heat was performed using the 
function C flt {T) = 5C an + C bg = ±5a™ + C bg , with 
a fit parameter 7 C , representing a generalized Griineisen 
parameter 26\. A parameterization of the form C hg — 
a + a\T + C12T 2 was used. The good quality of the fit, 
yielding SC an (cf. Fig.@| and a smooth background func- 
tion C hg (broken line in Fig. 13]), is reflected by the small 
residual (C(T) - C flt {T))/C(T) of less than 1.5 % in the 
range 4K < T < 10K. 

The notion of a second-order phase transition at 6 K is 
also corroborated by reassessing the low-T magnetic sus- 
ceptibility x(T)- The results of our experiment, being in 
good overall agreement with published data [H , are plot- 
ted in Fig. 2] as dx/dT vs.T. The measurements were 
performed by using a composite sample consisting of a 
large number of randomly oriented single crystals of total 
mass of 18.0 mg. As the figure indicates, dx/dT shows 
a distinct peak, the shape and position of which prac- 
tically coincide with the phase transition anomalies in 
oti(T) and C(T). Such a conformity between dx/dT and 



thermodynamic quantities, known from magnetic tran- 
sitions 31], indicates that spin degrees of freedom are 
involved to some extent in the transition. However, as 
will be shown below, their contribution alone appears in- 
sufficient to account for the observed entropy release. 

By integrating 5C an /T from 2 K up to 10 K, we obtain 
an entropy S an of (0.69 ± 0.05) J mor 1 K" 1 (cf. inset of 
Fig. [3]), with one mole substance containing 2-Na dimers 
(spins) and Na denoting Avogadro's constant. This rep- 
resents a lower bound of the true entropy S trans associ- 
ated with the phase transition, see below. The so-derived 
S an corresponds to about (6 ± 0.4)% of the system's 
full spin entropy (this would be i?ln2 for a hypothetical 
system containing Na spins with S = 1/2, with R the 
universal gas constant). It can be compared with calcu- 
lations of the magnetic entropy for 2D triangular-lattice 
5 = 1/2 systems, unfortunately available only for the 
isotropic (f = t) case. For the pure Heisenberg antifcr- 
romagnet [32j], which may serve as a rough estimate, we 
find a residual entropy for temperatures ksT < 0.04- J, 
corresponding to T < 10 K for J/ks = 250 K, of only 
(2.3 ± 0.1) % of i?ln2. By including ring exchange pro- 
cesses, likely to be significant near the Mott transition 
and which may account for the absence of Neel order- 
ing [HI , a somewhat l arge r value is expected. According 
to a variational study |12j | , these ring-exchange processes 
give rise to a low-T QSL state characterized by gapless 
fcrmionic excitations. This implies a spinon contribution 
to the specific heat of C spmon = j sp inon-T, with j spin on = 
(tt 2 /3)k%NAn(Ep) and a spinon density of states at the 
"Fermi surface" n(E F ) = 0.28/t spinon [12]. For t spinon 
= J [HI) we find ^spinon = 30 mJ mol _1 K -2 and a re- 
maining spinon entropy of 5.2 % of i?ln2 for T < 10 K. 
In Ref. 21] it has been argued that this QSL state is in- 
stable against the formation of spinon pairs - a nematic 
QSL state - by which S spmon is released to a full extent 
or partly (in the presence of impurity scattering). 

This residual spinon entropy S spmon appears to be 
quite close to the experimentally determined entropy re- 
lease S an , seemingly consistent with the scenario dis- 
cussed in ref. [2l[. However, by estimating possible er- 
rors, the discrepancy between S spmon and S an becomes 
significant. First, we note that the above-mentioned un- 
derestimation of S trans by S an is considerably large for 
the present case. The reason for this is mainly due to the 
peculiar non-mean-field shape of the transition, which 
required a truncation of the high- and low-T flanks in 
SC an (cf. Fig. HI). Simulations for quantifying the as- 
sociated error indicate that S trans / S an is probably in 
the range 1.2 - 1.4. Moreover, we expect that the avail- 
able low-T spin entropy will be further reduced for an 
anisotropic triangular lattice it' < t) - very likely the 
case for the present material 0, @|. Thus, an estimated 
entropy release S trans of about 7-8 % of Mn2 has to be 
contrasted with an available spin entropy of less than 5.2 
% of i?ln2. We consider this mismatch of about 30 - 50%, 
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FIG. 4: (Color online) Phase transition anomalies, derived as 
described in the text, in ct c and Qb (multiplied by a factor 
—2.5) (left scale), C (right scale), and d\/dT (scale is given 
by the arrow). 



by which the experimental value exceeds S spmon , to be 
a strong indication that spin degrees of freedom alone 
are insufficient to account for the phase transition. An 
obvious candidate to rationalize the mismatch would be 
charge degrees of freedom. In fact, according to optical 
conductivity measurements 33j, the charge-gap in this 
material is strongly suppressed indicating the material's 
close proximity to the insulator-metal transition. 

As a possible scenario we mention the proposed forma- 
tion of an exciton condensate . This is accompanied 
by a charge modulation breaking lattice symmetry. Ob- 
jections to this scenario may be due to the lack of clear 
signatures for a charge modulation below 6K in 13 C- 
NMR studies [16J. The poverty of evidence, however, 
might be due to the anomalous inhomogeneous broaden- 
ing of the NMR spectra 16], possibly obscuring a poten- 
tial line splitting. 

Alternatively, one might consider charge fluctuations 
at elevated temperatures, causing fluctuating electrical 
dipoles. Some evidence in support of this has been re- 
cently observed in dielectric measurements [34]]. Within 
such a scenario, the 6 K transition might indicate an or- 
dering of those preformed dipoles. We stress that a dis- 
tinct type of charge ordering, driven by an increase in the 
degree of frustration t' It , has recently been predicted for 
the present material 35}. Due to the ionic nature of the 
material, ordering phenomena in the charge sector are 
expected to give rise to distinct lattice effects (see, e.g. 
[23 . l3pl|). Since the charge distribution also affects the 
magnetic exchange constants via J oc t 2 /U with U the 
on-site Coulomb repulsion, a response in the magnetic 
susceptibility is also expected. 



In summary, measurements of the uniaxial ex- 
pansivities on the spin-liquid candidate k-(BEDT- 
TTF) 2 Cu2(CN) 3 yield distinct and strongly anisotropic 
lattice effects at GK. clearly identifying this phenomenon 
as a second-order phase transition. By means of a 
Griineisen-scaling Ansatz, the corresponding anomaly in 
the specific heat could be separated. Comparison of the 
associated entropy release with spin models suggests that 
charge degrees of freedom are involved in the transition. 
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partment of Energy Office of Science laboratory, operated 
under Contract No. DE-AC02-06CHII357. 



[1] 



[2] 



[3] 
[4] 
[5] 
[6] 
[7] 
[8] 
[9] 
[10] 

[11] 
[12] 
[13] 

[14] 
[15] 
[16] 
[17] 
[18] 
[19] 
[20] 
[21] 
[22] 
[23] 
[24] 
[25] 
[26] 
[27] 



[28] 
[29] 
[30] 
[31] 
[32] 

[33] 
[34] 
[35] 



N. Toyota, M. Lang, J. Miiller, Low- dimensional Molecu- 
lar Conductors in Springer Series in Solid-State Sciences 
154, Springer- Verlag Berlin Heidelberg 2007. 
A. Lebed (editor) The Physics of Organic Superconduc- 
tors and Conductors in Springer Series in Materials Sci- 
ence 110, Springer- Verlag Berlin Heidelberg 2008. 
Y. Shimizu et al, Phys. Rev. Lett. 91, 107001 (2003). 
W. Zheng et al, Phys. Rev. B 71, 134422 (2005). 
S. Yamashita et al, Nature Physics 4, 459 (2008). 
M. Yamashita et al, Nature Physics 5, 44 (2009). 
K. Nakamura et al, J. Phys. Soc. Jpn. 78, 083710 (2009). 
H.C. Kandpal et al, Phys. Rev. Lett. 103, 067004 (2009). 
T. Komatsu et al, J. Phys. Soc. Jpn. 65, 1340 (1996). 

G. Misguich et al, Phys. Rev. B 60, 1064 (1999). 

H. Morita et al, J. Phys. Soc. Jpn. 71, 2109 (2002). 
O.I. Motrunich, Phys. Rev. B 72, 045105 (2005). 

S. Yunoki, and S. Sorella, Phys. Rev. B 74, 014408 
(2006). 

Y. Hayashi et al, J. Phys. Soc. Jpn. 76, 053705 (2007). 
L.F. Tocchio et al, Phys. Rev. B 80, 064419 (2009). 
Y. Shimizu et al, Phys. Rev. B 73, 140407(R) (2006). 
Baskaran, Phys. Rev. Lett. 63, 2524 (1989). 
Kawamura et al, J. Phys. Soc. Jpn. 53, 4138 (1984). 
S. Lee et al, Phys. Rev. Lett. 98, 067006 (2007). 
Qi, and S. Sachdev , Phys. Rev . B 77, 165112 (2008). 
Grover et al, larXiv:0907.1710l 
Pott, and R. Schefzyk, J. Phys. E 16, 444 (1983). 
Miiller et al, Phys. Rev. B 65, 140509(R) (2002). 
Geiser et al, Inorg. Chem. 30, 2586 (1991). 
A. P. Ramirez et al, Physica B 280, 290 (2000). 

A. Briihl et al, Phys. Rev. Lett. 99, 057204 (2007). 
We exclude the possibility of a magnetic transition in- 
duced by sample inhomogeneities: no signatures for mag- 
netic order have been seen by NMR studies 0, [Hj]; no 
sample-to-sample variation has been found for the tran- 
sition temperature for the various crystals studied here. 
M. de Souza et al, Phys. Rev. Lett. 99, 037003 (2007). 
M. de Souza et al, Phys. Rev. Lett. 101, 216403 (2008). 
M. de Souza et al, Physica B 404, 494 (2009). 

M.E. Fisher et al, Phil. Mag. 7, 1731 (1962). 

B. Bernu and G. Misguich, Phys. Rev. B 63, 134409 
(2001). 

I. Kezsmarki et al, Phys. Rev. B 74, 20110RR) (2006). 
M. Abdel-Jawad et al, unpublished results. 

H. Li et aUarXiv:0908.4109 



